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Abstract. We present full-star simulations of Type la supernova explosions on the basis of the standard Chandrasekhar-mass 
deflagration model. Most simulations so far considered only one spatial octant and assumed mirror symmetry to the other 
octants. Two full-star models are evolved to homologous expansion and compared with previous single-octant simulations. 
Therefrom we analyze the effect of abolishing the artificial symmetry constraint on the evolution of the flame surface. It turns 
out that the development of asymmetries depends on the chosen initial flame configuration. Such asymmetries of the explosion 
process could possibly contribute to the observed polarization of some Type la supernova spectra. 
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1. Introduction 

Recent observations have shown that spectra of — at least 
some — Type la supernovae (SNe la) show significant polariza- 
tion (e.g. Ho well et all l2001: Wa ng et all l2003). The origin of 
this effect is stil l unclear, but some possible mechanisms hav e 
been proposed JWang et all l2003t iKasen et all l2003l |2004) . 
Besides peculiar spatial structures of the remnant and interac- 
tion of the ejecta with the companion star, anisotropics of the 
explosion process itself could contribute to the observed polar- 
ization. 

This idea fits well into the currently favored astrophysi- 
cal model of SNe la. These events are assigned to thermonu- 
clear explosions of carbon/oxygen white dwarf (WD) stars 
in binar y systems (for a review of S N la explosion mod- 
els see Hilleb randt & Niemevetl l2000l) . From non-degenerate 
companions they accrete matter until they approach the 
Chandrasekhar mass. After some thousand years of convective 
carbon burning thermonuclear reactions in the center of the star 
form a flame, i.e. a small spatial region undergoes thermonu- 
clear runaway. This flame can travel outward in two modes. 
If it is mediated by shock waves, the propagation proceeds at 
supersonic velocities leading to so-called detonations, while a 
mediation due to heat conduction of the degenerate electrons 
leads to subsonic deflagrations. A prompt detonation can be 
excluded as a model for SNe la si nce it fails t o produce the 
observed composition of the ashes ( Arnett, 1969). It burns the 
entire star at high densities producing mainly iron group ele- 
ments. 
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Consequently a mechanism has to be developed wherein 
the star can expand prior to incineration. This is provided if 
the flame propagation starts out as a deflagration. This flame 
propagation mode, however, is much to slow to explain SN la 
explosions. Realizing that the astrophysical scenario is char- 
acterized by Reynolds numbers as high as 10 14 , one expects 
the development of strong turbulent features in the flow. These 
will interact with the flame wrinkling and stretching its sur- 
face. Since the fuel consumption rate of deflagration flames is 
determined by the flame surface area, turbulence will signifi- 
cantly accelerat e the flame pr opagation. Three-dimensional SN 
la models (e.g. Hillebrandt et al., 2000) taking this effect into 
account have been shown to lead to ro bust explosions and re- 
produce many featu res of SNe la (e.g. iReinecke et al1 Eo02b: 
Ga mezo et all2003l) . There have been speculations on the pos- 
sibility of a later transition from the deflagration to a detona- 
tion, but a viable mechanism to trigger this transition in SNe la 
explosions could not be identified yet In any 

case, the deflagration stage requires thorough investigation es- 
pecially in the context of asymmetries of the explosion process 
for it is intrinsically anisotropic. Due to buoyancy (Rayleigh- 
Taylor) instabilities the large scale features of the flame show 
rising "bubbles" of hot and light ashes and in between down- 
drafts of cold and dense fuel. It is still an open question, how 
strong the resulting asymmetries may become. This issue will 
be addressed in the present study. 

It is clear that in order to explore asymmetries of the explo- 
sion process it is mandatory to simulate the full star. Most of 
the three-dimensional simulations so far have been restricted to 
one spatial octant assuming mirror symmetry to the other oc- 
tants. Obviously, this approach imposes an artificial symmetry 
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constraint on the flame evolution and our goal in this work is 
to study the effects that result from removing this symmetry. 
There has been only one previous simulation of the deflagra- 
tion phase in the full star dCalderet all 12004). but it failed to 
explode. The reason for this is presumably that it was started 
with very artificial conditions for the initial flame (a perfect 
sphere). Nevertheless, the merit of this model is that it points 
to the importance of the initial flame configuration in the con- 
text of asymmetries of the explosion process. 

The flame ignition is difficult to address in numerical sim- 
ulations and afflicted with large uncertainties (iWooslev et all 
120041) . We will explore the impact of initial flame configura- 
tions on the explosion models in a parametrized way. As a 
first step, in this study we investigate anisotropies in the flame 
evolution resulting from simulations covering the full star in- 
stead of only one octant and from moderately asymmetric ini- 
tial flame configurations. This allows comparison with previous 
simulations. Detailed studies of multi-spot ignition scenarios 
and strongly dipolar ignited models are subject to forthcoming 
publications. 

2. Implementation 

2.1. Flame Model 

The stratification of the hot and light ashes below cold 
and dense fuel leads to a buoyancy instability known as 
the Rayleigh-Taylor instability. In its non-linear stage, bub- 
bles of burning materi al rise into the fuel with a velocity 
JPavies & TavlorHl950l) 

V oc y/gLAt, (1) 

where g and L denote the gravitational acceleration and the 
length scale of the bubble, respectively. The Atwood number 
At characterizes the density contrast between the interior and 
the exterior of the bubbles. At the interfaces of the burning bub- 
bles shear flows give rise to a secondary instability, the Kelvin- 
Helmholtz instability. It generates turbulent velocity fluctua- 
tions at this (integral) length scale L, which decay to smaller 
scales in a turbulent cascade over a so-called inertial range and 
are dissipated at the Kolmogorov length. For the velocity fluc- 
tuations v at some scale I in the inertial range Kolmogorov- 
scaling yields 

/A 1/3 

v(t) = v(L)l-\ . (2) 

The resulting turbulent motions interact with the flame defining 
the deflagrations in thermonuclear supernova explosions as a 
problem of turbulent combustion. 

The correct modeling of turbulent combustion seems to be 
the crucial issue in deflagration models of SNe la. These can 
greatly benefit from the substantial progress that has been made 
in turbulent combustion research over the past decades. A fun- 
damental notion is that turbulent combustio n proceeds i n differ- 
ent regimes (see Damkohler ( 1940) and Borghi ( 1985); Peters 
( 1986. l2000l) for more recent work). The so-called flamelet 
regime of turbulent combustion applies to most parts of the de- 
flagration in thermonuclear supernovae. 



In this regime, turbulent eddies wrinkle and stretch the 
flame but do not penetrate its internal structure. The internal 
flame structure (which extends to less than a centimeter) as well 
as the smallest scales of flame interaction with turbulence are 
tiny compared with the dimensions of a WD star (~10 8 cm). It 
is therefore well justified to treat the resulting "turbulent flame 
brush" as a discontinuity (the so-called flame front) separating 
the fuel from the ashes, as long as the small-scale effects are 
appropriately modeled. 

Two approaches are taken to ensure that. First, known ef- 
fects have to be taken into account. This applies mainly to the 
interaction of the flame with turbulence on unresolved scales. 
Wrinkling will increase its surface and therefore enhance the 
burning rate. The shortcoming in resolving the flame surface 
can be compensated by assigning an increased turbulent burn- 
ing velocity to the flame front reproducing the correct fuel 
consumption rate. It turns out that in the flamelet regime the 
turbulent burning velocity is proportional to the turbulent ve- 
locity fluctuations and decouples completely from the laminar 
burning speed determined by the microphysics. Second, one 
may suspect that on unresolved scales n ew physical effects oc- 
cur (l ike "active turbulent combustion". iNiemever & WooslevI 
1 997 ) that affect the flame propagation on large scales. This has 
to be studied in separate small-scale models focusing on a nar- 
row win dow in scale space. Simulations of the cellular burning 
regim e jRopke et all 120031 IRopke et all l2004albt bell et all 
2004) led to the conclusion that the resulting increased burning 
velocity can be "renormalized" with little effect on the large 
scale flame propagation in most stages of the SN la explosion. 
At least from this regime drastic effects are unlikely. 

Due to the expansion of the WD, the density at which burn- 
ing takes place steadily decreases. This is accompanied by a 
gradual broadening of flame width. At some point turbulent 
eddies can affect the flame's internal structure and burning en- 
ters the so-called distributed regime. Since this is expected to 
happen at very late stages, it may not contribute much to the 
production of iron group nuclei. Nevertheless, it could signifi- 
cantly increase the abundances of intermediate mass elements 
and hence the total energy release with consequences for mod- 
eling of lightcurves and spectra. This regime is not yet accessi- 
ble in our code but current efforts aim at including it in future 
simulations. In this respect the results presented here are pre- 
liminary. 

To summarize the turbulent flame propagation in thermonu- 
clear supernova explosions we note that large-scale features are 
driven by buoyancy according to Eq. This sets the integral 
scale of a turbulent cascade to smaller scales where the velocity 
fluctuations in the simplifying picture of Kolmogorov scaling 
follow Eq. 0. These velocity fluctuations determine the turbu- 
lent flame velocity which thus falls off slower toward smaller 
scales than the square root law of Eq. {Q. 

2.2. Numerical methods 

The implementation of our astrophysica l model in a numerical 
scheme uses the techniques described bv lReinecke eTaD {l 999 , 
2002a). Only a brief outline will be given here. 
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Fig. 1. Initial flame shape of model c3_4n. 



The hydrodynamics is modeled by the reactive Euler 
equations. Their numerical treatment is based on the 
Prometheus implementatio n dFrvxell et all Il989l) of the 
Piecewise Parabolic Method JColella & Woodwardlll984f) and 
augmented by an equation of state for WD matter. The descrip- 
tion of th e nuclear reactions follow s the simplified scheme pro- 
posed by iReinecke et al.l J2002ab and comprises five species: 
a-particles, 12 C, ls O, "Mg" as a representative of intermediate 
mass elements, and "Ni" representing the iron group. 

The propagation of the turbulent deflagration flame ismod- 
eled applying the level set technique dOsher&Sethiani ri988). 
Its basic idea is to associate moving fronts with the zero level 
set of a signed distance function G, which is then propagated 
in an appropriate way. A modified version of this technique has 
been successfully used to model terrestrial deflagration flames 
Smi lianovski et alJ dl997l) . We apply the so-cal led passive im- 
pleme ntation of the method (for details see IReinecke et all 
1999). The propagation of the flame front via level sets re- 
quires knowledge of the flame propagation velocity, which is 
composed of two parts: the advection of the flame in the flow 
and its self-propagation due to burning. The former is handled 
by the hydrodynamics scheme while for the latter the G-field is 
updated according to the turbulent burning velocity. To deter- 
mine this quantity we apply a sub-grid scale turbulence mo del 
first a dopted to SN la simulations by Niemever & Hillebrandt 
(19951). 



Our approach is similar to Large Eddy Simulations (LES) 
known from computational fluid dynamics. The efficiency of 
the numerical scheme is a result of the combination of level- 
set techniques with a sub-grid scale model. It enables nu- 
merically converged simulations at ra ther crude resolutions 
( IReinecke et alll2002ald:fRopkeL Eo04). This makes parameter 
studies based on thre e-dimensional simulations possible (e.g. 
Ropke & HillebrandiHoQl. 



Besides the wide range of spatial scales, the significant ex- 
pansion of the WD during the explosion complicates numeri- 
cal simulations of thermonuc lear supernovae. To ca pture this 
expansion, the simulations by Reinec ke et all J2002albld) were 
carried out on a nonuniform Cartesian computational grid with 
fine resolution in the inner parts and an exponentially expand- 
ing grid spacing in the outer parts of the domain. This, however, 
causes very low resolution of the flame front in late stages, es- 
pecially i n direc t ions d iagonal to the axes. To overcome this 
problem, Ropke (2004) developed a scheme to simulate ther- 
monuclear supernovae on an non-static uniform Cartesian com- 
putational grid which tracks the expansion of the WD. Besides 
a number of numerical advantages it offers a possibility to 
follow the evolution of the explosion to homologous expan- 
sion with only moderate additional computational expenses. 
Furthermore, at a given time, it provides a constant resolu- 
tion throughout the domain. This is important in full-star sim- 
ulations without artificial symmetries. On a non-uniform grid 
main features of the flame front could easily propagate in re- 
gions with poor spatial resolution. This is avoided in the co- 
expanding uniform grid approach. For detailed description and 
thorough testing of this scheme we r efer to lRopkel (Eo04). 



2.3. Simulation setup 

The setup used in the simulations presented below is similar to 
that described bv lRonkel 12004). except for the computational 
domain now covering the full WD star instead of only one oc- 
tant. All simulations were carried out on a co-expanding uni- 
form Cartesian computational grid of [512] 3 cells. The initial 
resolution was chosen to be Ax = 7.9 km. Reflecting bound- 
aries were applied on all sides of the domain. 

The simulations started with a cold isothermal WD star of 
1.406 M Q with a central density of p c = 2.9 x 10 9 gcirT 3 and 
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Fig. 2. Initial flame shape of model fl. In the projections, the solid lines indicate the center of the WD, while the dashed lines 
indicate the center of the flame configuration (the ^-coordinates coincide at the chosen scale of the plot). 



a temperature of To — 5 x 10 5 K. Its initial composition was 
chosen as a 50% mixture of carbon and oxygen. In this config- 
uration we placed the initial flame by instantaneously burning 
the material behind the front. 

2.4. Initial flame configuration 

A crucial parameter of thermonuclear supernova simulations is 
the prescription of the flame ignition. This problem is, how- 
ever, difficult to tackle in explosion simulations. The formation 
of the initial flame front(s) is determined by about a thousand 
years of convective burning prior to thermonuclear runaway. 



There have be en only few attempts to investigate the igni- 
tion c onditions JGarcia-Senz & Wooslev[ll995llWooslev et all 
2004) and these conclude that flame formation is likely to pro- 
ceed in a number of burning "bubbles" that may rise to dis- 
tances up to 200 km from the center before the runaway. 

Initial flame configurations have to be represented on the 
rather coarse computational grids used in multidimensional 
thermonuclear supernova simulations. Several possibilities for 
this have been suggested in literature. In the simplest approach 
one start s out with a fl ame front in form of a single sphere 
(e.g. iGamezo et all 120031) . In this case, however, the high ini- 
tial symmetry is broken by numerical noise and it takes a 
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substantial amount of computer time until the first Rayleigh- 
Taylor features develop. This is probably not a very realistic 
scenario, since the multitude of burning bubbles predicted by 
Woosl ev et al 1 ll2004 provides seeds for the development of 
nonlinear buoyancy instabilities from the very beginning of the 
flame propagation. Therefore one may start out with a sphere 
on which perturbation s are imposed, such as the c3 model of 
Reinecke et al .1 (12002 a). A better approximation is provided by 
multi-spot ignition scenarios IReinecke et all l2002bh . where 
some burning bubbles are initialized which, due to the lack of 
resolution, represent a collection of real flame spots. Here, the 
problem is, however, that a rather high resolution is required to 
accommodate a decent number of burning bubbles (note that 
the high resolution is not necessarily required for numerical 
convergence). The models in this paper comprise [256] 3 com- 
putational cells per octant rendering reasonable multi-spot igni- 
tion scenarios impossible. Therefore we chose to represent the 
initial flames in one of our simulations by a cluster of bubbles. 
These overlap in the center while some of them float detached 
from the bulk. The result is a clumpy or foam-like initial flame 
structure which we think provides a reasonable approach given 
the restricted computational resources. 

Since the initial flame configurations turn out to be a crucial 
parameter for the simulations, we will now describe our models 
in detail. 

As a standard test case IReinecke et al J J2002al) established 
the c3 model which is very simple and allowed t he comparison 
of two - and three-dimensional simulations (see lReinecke et all 
2002a). Here the initial flame is represented by a sphere which 
is perturbed by a superposition with toroidal rings. This config- 
uration is visualized in Fig.^ where the isosurface corresponds 
to G = representing the flame front. We will denote the full- 
star version presented here c3_4n. Obviously, the center of the 
initial flame coincides with the center of the WD in this model. 

A second initial flame configuration was obtained by the 
superposition of a number of bubbles around the center of the 
WD star. The radius of the bubbles was set to 3.5 km and their 
locations were chosen in a random selection process under the 
constraint of a Gaussian distribution around the center with 
a dispersion of <x = 100 km. The resulting configuration is 
depicted in Fig. |2] We will refer to the corresponding model 
as fl in the following. Note that the flame distribution result- 
ing from the random selection of locations is not symmetric 
and that the outermost flame spots are located at a distance 
of ~20 km from the center in agreement with IWooslev et alJ 
Resulting from the random choice of a limited num- 
ber of "bubbles", the center of the initial flame is displaced by 
r = (0.3, -13.0, 3.3) km from the center of the WD, i.e. it is 
located at a distance of 1 3 .4 km off-center. 



3. Morphology of the flame evolution 

Analyzing the results of our simulations, we will first describe 
the flame evolution in our models. Global characteristics of the 
explosions will be discussed in Sect.|4] 




Fig. 3. Comparison of the c3_4n full star simulation (light gray) 
with the c3_oct single-octant simulation mirrored to the other 
octants (dark gray) at t = 1 .5 s with an edge length of the cube 
of 2. 1 6 x 1 9 cm (top) and at t — 10.0s with and edge length of 
the cube of 2.40 x 10 10 cm (bottom) 



3.1. Effects from removing the symmetry constraints 

Ropkd d2004 applied a c3 flame ignition condition for a 
demonstration of the co-expanding grid approach in three di- 
mensions. There only one octant was simulated assuming mir- 
ror symmetry to the other octants. We will now compare this 
simulation (denoted c3_oct) with our c3_4n full-star simula- 
tion. Although the c3 initial flame shape is not a realistic sce- 
nario, this model offers a way to assess possible differences in 
the full star implementation. 

A comparison of the flame shapes at t = 1 .5 s and at 
t = 10.0 s is given in Fig. [3] The isosurfaces again represent 
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G — 0, which is associated to the flame front at early times, 
while at later times when the burning has ceased it still pro- 
vides a g ood approxima tion to the interface between fuel and 
ashes (cf. Ropke, 2004). Fig. [3] shows that the evolutions of 
the flame structures of the c3_oct simulation (here mirrored to 
the other octants) and the full-star model proceed very simi- 
lar. The G = isosurface in the full-star simulation is in close 
agreement with the octant symmetry of the c3_oct model, even 
after an evolution time of 10 s. Obviously, the axial symme- 
try of the initial flame front is broken in the same way in the 
two implementations. It turns out, that this happens preferen- 
tially at locations where the flame front is aligned diagonal to 
the computational cells of the Cartesian grid indicating that the 
symmetry breaking is induced by discretization effects. 

It seems appropriate to point out that the deflagration SN la 
model does not develop new large-scale features when the com- 
putational domain is extended to the full star. Unlike convec- 
tive burning pheno mena (e.g. pre-ignitio n evolution of a WD 
toward a SN la, see Wooslev et al. 2004, or core-collapse su- 
pernovae, see Schec k et al.ll2004l) our models form no global 
dipole-asymmetries. The reason for this is that the rapid prop- 
agation of the burning front does not allow several eddy 
turnovers which seem to be crucial for the development of such 



configurations. We also find no indication for the development 
of larger plume structures in the c3_4n simulation than in the 
c3_oct model. 

Thus, as long as the initial conditions (especially the initial 
flame configurations) do not violate the mirror symmetry, it is 
safe to restrict the simulations to a single spatial octant. Also 
for our full-star models the result is encouraging since it proves 
that they are able to preserve initial symmetries, without explic- 
itly imposing those by symmetries of the domain. This verifies 
our full-star implementation. Hence asymmetries as described 
below are of physical nature and not a numerical artifact. 

3.2. Effects of asymmetric initial flame configurations 

We will now address the question how the SN la explosion pro- 
ceeds if the initial flame configuration does not possess any 
symmetry. This is the case for our/i model. Starting from the 
initial flame configuration shown in Fig. [2] the evolution of the 
flame front in the explosion process is illustrated by snapshots 
of the G = isosurface at t — 0.3 s, t = 0.6 s, and t = 2.0 s in 

Kg. a 

The development of the flame shape from ignition to 
t = 0.3 s is characterized by the well-known "mushroom-like" 
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Fig. 6. Model/7 at t — 10 s. In the projections the column mass density is linearly color-coded. 



structures resulting from buoyancy. This is especially well vis- 
ible for the bubbles that were detached from the bulk of the 
initial flame. But also the perturbed parts of the contiguous 
flame closer to the center develop nonlinear Rayleigh-Taylor 
like features. This flame evolution generates strong vortices in 
the resolved flow field, illustrated by the temporal evolution 
of the velocities in a two-dimensional slice in Fig. [5] The or- 
der of magnitude of the resulting velocities is in agreement 
with Eq. Q, given the length scale of the buoyant structures 
of ~100km, a gravitational acceleration of about 10 I0 cms~ 2 , 
and an Atwood number of the order of 0. 1 . Due to the multitude 
of ascending structures, the flow pattern is strongly irregular. 



During the following flame evolution inner structures of 
smaller scales catch up with the outer "mushrooms" and the 
separated structures merge forming a more closed configura- 
tion (see snapshot at t — 0.6 s of Fig.0}. This is a result of the 
large-scale flame advection in the turbulent flow and the expan- 
sion of the ashes. 

Up to this stage the flame was strongly anisotropic. 
However, in the later evolution a preferentially lateral growth of 
bubbles filled with ashes smooths out parts of the anisotropies. 
As this effect continues to act even after the burning has ceased 
in our model, it is of hydrodynamic origin rather than caused by 
flame propagation. The flame develops a more spherical shape 
and only a slight anisotropy is retained. This is apparent from 
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Fig. 7. Temporal evolution of the total energies in our models. 



the snapshot at t — 2.0 s in Fig.0] At the stage depicted there 
self-propagation of the flame due to burning has terminated. 

The subsequent evolution is characterized by the approach 
to homologous expa nsion which i s reached to a good approx- 
imation at t ~ 10 s (Rtinke, 2004). A snapshot corresponding 
to t — 10 s is shown in Fig.|6] The differences to the snapshot 
at t — 2.0 s in Fig. [4] are small, but the smoothing effect has 
rounded the structure somewhat further. This confirms the as- 
sumption that it is caused by the hydrodynamics. 

It is interesting to note that the fl model, though it possesses 
a slight large-scale asymmetry, finally develops a smoother sur- 
face of ashes than the c3_4n model (compare the snapshots at 
t - 10 s in Figs.[3]and[6j. This can be attributed to the fact that 
the amplitude of large-scale perturbations (interpreted as devi- 
ations from a spherical shape) in the fl initial flame is smaller, 
although it is more perturbed on small scales. 

An important aspect of full-star simulations are the result- 
ing asymmetries in the density distribution. An impression of 
these is given by the projections of the density to the coordinate 
planes in Fig. [6] The central high-density region is aspherical 
and also in the outer regions filaments of increased density are 
visible. These result from downdrafts of dense fuel between 
bubbles of light ashes during the explosion. 

Despite the slight initial misalignment of the centers of the 
WD and the flame, we observe no buoyancy-driven ascension 
of the flame as a whole toward the surface. For this a large-scale 
flow around the flame structure would be necessary. Such a 
flow-structure can, however, not develop since our initial flame 
configuration causes an irregular vortical velocity field, which 
dominates the flame advection. We emphasize here, that this 
large-scale effect is not caused by the self-propagation of the 
flame and should thus be independent of the particular turbu- 
lent flame model. 



Fig. 8. Temporal evolution of the turbulent energies in our mod- 
els. 



Table 1. Global characteristics of the models. 



model 


m ("Ni") 


m("Mg") 


Etot 


P 

"nuc 




[M ] 


[M Q ] 


[10 50 erg] 


[10 50 erg] 


c3_4n 


0.507 


0.151 


4.275 


9.466 


fl 


0.579 


0.165 


5.585 


10.776 



4. Global parameters of the explosion models 

The global characteristics of the models are summarized in 
Tabled It is obvious that the/7 initial flame configuration leads 
to a more vigorous explosion than the c3_4n model. It produces 
14% more iron group elements and consequently the release of 
nuclear energy is 14% higher. This amounts to a 31% increase 
in the total explosion energy. The excess in produced interme- 
diate mass nuclei is only 9% and thus the relative fraction of 
these is lower in the fl model. The temporal evolutions of the 
total energies in our simulations are plotted in Fig. [7] 

Two aspects of the fl model contribute to the more ener- 
getic explosion. First, the initial flame configuration extends 
to slightly larger radii than the c3_4n initial flame (compare 
Figs. Hand GJ Therefore the outer flame fronts experience a 
larger gravitational acceleration and the development of struc- 
tures resulting from buoyancy instabilities is accelerated. This 
leads to an enhanced generation of turbulence and conse- 
quently the value of the turbulent flame propagation velocity 
is higher. Hence more material is burnt at higher densities re- 
sulting in iron group nuclei. 

Second, due to its composition of small bubbles the ini- 
tial flame shape of the fl model is perturbed with shorter 
wavelengths. This provides more seeds for the development of 
buoyancy instabilities and consequently the wrinkling of the 
flame is is enhanced resulting in an increased fuel consump- 
tion. Contrary to that, in the c3_4n model it takes some time 
until the axial symmetry is broken and "mushroom-like" struc- 
tures form. The fl model therefore seems more appropriate to 
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Fig. 9. Angular averaged distribution of the ejecta of model Fig. 10. Angular averaged distribution of the ejecta of model/7 
c3_4n in radial velocity space after 10 s. in radial velocity space after 10 s. 



mimic realistic initial flame structures (perturbed by convective 
motions and eventually igniting in multiple spots) on coarse 
computational grids. The evolutions of the turbulent energies 
in both simulations are plotted in Fig. [S] showing clearly the 
enhanced turbulence production in the fl model. Since the tur- 
bulent burning velocity scales with the turbulent velocity fluc- 
tuations, the flame propagates faster in that model. 

Synthetic light curves and spectra derived from explosion 
models are largely determined by the distribution of the ejecta 
in velocity space. However, the interpretation of these data 
is not straightforward and requires caution. One of the for- 
mer main uncertainties in analyzing the results from three- 
dimensional explosion simulations, namely the short simula- 
tion time, could be overcome in the models presented here. 
On the co-expanding computational grid they were evolved 
to 10 s providing a fair approximation to homologous expan- 
sion (Ropke, 2004). However, most of the currently available 
numerical codes treat the radiation transport one-dimensional. 
One should keep in mind, that multidimensional structures 
(cf. FiglfJ could considerably change the results from three- 
dimensional synthetic light curve and spectra calculations de- 
pending on the line of sight. 

In Figs.l9landll0lwe present the angular averaged profiles 
of the velocity distribution of the ejecta. These figures show 
clearly that both explosion models are too weak to accelerate 
the outer ejecta to the highest observed velocities (see for ex- 
ample the velocitie s of the Si n, S n and Ca n lines given by 
iBenetti etaill2004 . While they reach to -1.2 x 10 9 cms _1 in 
the c3_4n model, the outermost ejecta of the fl model are found 
at velocities of 1.5 x 10 9 cms _1 . 

Due to the close similarity of the c3_4n simulation to the 
c3 model calculated in only one octant, the distribution of the 
species in the ejecta shows no visible differenc es (com pare with 
the analogous plot for the c3 model of Ropke, 2004). However 
the dev iations from the the profiles given by Rei necke et alJ 
(2002b) are significant, especially at higher velocities. This is 
simply due to the fact, that we evolved the model to 10 s, while 
former simulations stopped at 1.5 s where homologous expan- 



sion is not yet reached (Ropke, 2004). At t - 1.5 s the distri- 
bution of the elements is quite uniform, but this is not retained 
in the approach to homologous expansion. The later accumu- 
lation at lower velocities can probably be attributed to the still 
non-negligible gravitational potential. 

A feature of the c3 model seems to be that little fuel is left 
at the center of the star. For this, a possible explanation would 
be that with the given initial flame configuration it takes longer 
to break the axial symmetry. True three-dimensional features 
with "mushroom-like" bubbles of ashes and downdrafts of fuel 
form later. Therefore less unburnt material sinks down to the 
center. 

This is different in the species distribution of the fl model 
shown in Fig.[l0] Here the higher amount of unprocessed ma- 
terial near the center can be attributed to the accelerated devel- 
opment of the typical buoyancy-induced flame structures. The 
earlier formation of downdrafts transports unburnt matter more 
efficiently toward the center. However, this effect may be partly 
compensated when the model is started fro m a true multi- 
spot i gnition scenario (cf. to the b30 model of iTravaelio et all 
2004]). In this case the increased surface of the flame near the 
center leads to a more efficient burning. A possible temporal 
distribution of the formation of initial flames would also favor 
an enhanced burning near the center. 

A second difference between our two models is the distri- 
bution of iron group nuclei (denoted as "Ni"). The accelerated 
flame evolution and the stronger explosion mix these elements 
further out in the fl simulation. A similar, but weaker effect 
can be observed for the intermediate mass nuclei ("Mg"). The 
low production of intermediate mass elements in both mod- 
els is in part due to the small computational grid applied in 
the simulations. Naturally, subsequent stages of the burning are 
less resolved in t he co-expand ing grid approach. By means of a 
resolution study Ropke (2004) found that starting with [256] 3 
cells per octant provides sufficient resolution in the early stages 
of flame propagation when iron group nuclei are synthesized. 
However, the later production of intermediate mass elements 
may be underestimated by about 50% due to the coarse reso- 
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lution of the flame. This, in turn, leads to a somewhat lower 
total energy of the models. Also the lower relative production 
of intermediate mass elements in the fl simulation may par- 
tially be caused by a resolution effect. Since this model ex- 
plodes more vigorously, the WD expands faster. The computa- 
tional grid tracking the expansion will therefore be coarser at 
the stages of intermediate mass element production in the fl 
model than in the less energetic c3_4n simulation. 

5. Conclusions 

We presented two full-star simulations of deflagration models 
for thermonuclear supernova explosions. These were started 
with different initial flame configurations. Both simulations re- 
sulted in robust explosions, though they were rather on the 
weak side compared with observations. The models allowed 
an analysis of the effects of removing the artificial symme- 
tries of former simulations. The c3_4n model retained the ax- 
ial symmetry of our standard c3 example although it was cal- 
culated on a full-star domain. We note that the abolishment 
of the domain-symmetries alone does not introduce any addi- 
tional asymmetries to the model. This result validates previ- 
ous single-octant simulation s (e.g. iReinecke et all l2002albTcl 
Ron ke & Hiilebrandl 120041) since it shows that the symmetry 
constraint imposed there did not suppress physical effects. 

To introduce asymmetries into our models it is necessary to 
start with an asymmetric initial flame configuration. This was 
realized by the "foamy" initial flame surface of the fl simula- 
tion. Starting the flame propagation with such an anisotropic 
configuration leads to asymmetries in the ejecta. However, the 
initial asymmetries are partially smoothed out by the expan- 
sion of the ash-filled bubbles, which in late phases proceeds 
preferentially in lateral directions. To what degree the retained 
asymmetry could contribute to the polarization of the spectrum 
remains to be tested by multidimensional radiative transfer cal- 
culations. 

As a consequence of the perturbed initial flame shape an 
irregular flow pattern establishes. This has strong impact on 
the flame evolution. In the/7 simulation two effects are noted. 
First, small bubbles in the initial flame configuration that were 
detached from the bulk merge with it shortly after ignition re- 
sulting in a more or less connected flame structure. Second, the 
model is robust against slight misalignments of the centers of 
the flame and the WD. No buoyancy driven rise of the flame 
as a whole toward the sur face of the star wa s observed. This is 
in contrast to the result of lCalder et al.l 1120041) . who found that 
the flame quickly ascends from the center although the initial 
misalignment was somewhat smaller t han in our fl model. The 
reason thereof seems to be that Calder et al. (2004) started with 
flame in shape of a perfect sphere. This configuration possesses 
no seeds for the development of Rayleigh-Taylor instabilities 
which consequently have to be introduced by numerical effects. 
Therefore the establishment of a vortical velocity field is de- 
layed and a symmetric large-scale flow may develop driving the 
flame to larger radii. Differences in the flame models could also 
effect the results, but seem unlikely t o be th e primary reason for 
th e discrepancy. W hile Calder et al. (2004) use the flame model 
of lGamezo et alJ J2003I) which is based on the assumption of a 



buoyancy dominated scaling according to Eq. {!}, we model 
the turbulent flame speed taking into account the interaction of 
the flame with turbulence on unresolved scales. From Eq. it 
is clear that this results in a slower decrease of the flame speed 
toward smaller scales. 

In this first study only one anisotropic initial flame config- 
uration was tested. Open questions remain regarding the flame 
evolution in case of stronger initial asymmetries, especially 
more pronounced dipole shapes. How robust are the models 
to stronger displacements of the initial flame from the center? 
What asymmetries can then result from deflagration models? 
An extre me case would be a on e-sided ignition which was sug- 
gested bv lWooslev et alJj2004h . Such models will be addressed 
in a forthcoming study. 
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